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ABSTRACT 

Damping mechanics for simulating the damped dynamic characteristics in specialty 
composite structures with compliant interlaminar damping layers are presented. Finite- 
element based mechanics incorporating a discrete layer (or layer-wise) laminate damping 
theory are utilized to represent general laminate configurations in terms of lay-up and fiber 
orientation angles, cross-sectional thickness, shape and boundary conditions. Evaluations of 
the method with exact solutions and experimental data illustrate the accuracy of the 
method. Additional applications investigate the potential for significant damping 
enhancement in angle-ply composite laminates with cocured interlaminar damping layers. 


*NASA Resident Research Associate at Lewis Research Center. 



INTRODUCTION 


The passive damping capacity of polymer matrix composites is receiving current attention 
[1]. Composite materials are readily preferred in many flexible structures, where passive 
damping may significantly improve dynamic and acoustic performance. Previous research 
has shown, that composite materials can provide low but critical levels of damping 
[2,3,4]. Yet, due the inhomogeneity in properties and anisotropy through the thickness 
of laminated composites, additional enhancements of passive damping are possible from 
contributions of interlaminar shear damping in specialty composite laminates with 
embedded interlaminar damping layers [5,6,7,8,9,10]. 

Most of the previously mentioned work has provided valuable insight to the damping of 
specialty composite laminate with damping layers, yet, it involves limitations which restrict 
applications to composite structures of general lay-up, geometry and boundary conditions. 
Consequently, this paper describes structural damping mechanics which enable the 
prediction of passive damping and other dynamic characteristics, such as, natural 
frequencies, mode shapes and transfer functions in specialty composite structures with 
embedded damping layers. The mechanics are generalized to include effects of general 
ply orientations, lay-ups, cross-sectional thickness, structural configurations and boundary 
conditions. A discrete layer damping theory (DLDT) is utilized with the capability to 
represent damping contributions from extension, flexure, flexure-shear coupling, and 
coupling from other possible asymmetries that may exist due to the heterogeneous nature 
of the laminate. A novel finite element based method incorporating DLDT is described 
and utilized to model the damped characteristics of composite structures of general 
thickness, shape, lamination, and boundary conditions. In this context, a four-node 
bilinear plate finite element is developed and described herein. 
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Damping and natural frequency predictions are compared with both exact solutions for 
simply supported cross-ply composite plates [10], as well as, with experimental data for 
unsupported angle-ply plates to illustrate the quality and versatility of the damping 
mechanics. Additional evaluations and applications investigate the potential to enhance 
the material damping in angle-ply laminates and the associated trade-offs in stiffness and 
natural frequencies. Some of the results may be also used to establish guidelines for 
damping measurement and damping improvement in these specialty composite structures. 
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MECHANICS FOR COMPOSITE DAMPING 


The current section overviews the composite damping mechanics enabling the analysis of 
composite structures engineered for high passive damping. 

Composite Plies 

For a unidirectional composite loaded along the material axes, micromechanical models 
have been developed for the synthesis of elastic [11] and dissipative properties [12]. 

Five independent elastic parameters completely characterize the stiffness of a 
unidirectional composite (orthotropic but transversely isotropic material): in-plane 
longitudinal modulus E, n , transverse in-plane modulus in-plane Poisson’s ratio v, 12 , 
in-plane shear modulus G 112 , and interlaminar shear modulus G^. Additionally, four 
independent damping loss factors characterize the composite damping, that is, 
longitudinal loss factor rj n (direction 11), transverse in-plane damping rj a (direction 22), 
transverse through-the-thickness damping ijb-Ub (direction 33 ), in-plane shear damping 
r Jl6 (direction 12 ), interlaminar shear damping r? 14 (direction 23 ), and interlaminar shear 
damping 17*= r? I6 (direction 13 ). Both damping and elastic properties are explicitly related 
to fiber/matrix properties, and the fiber volume ratio (FVR). 


For the case of off-axis composites, ie. composites loaded at an angle 0 , the stiffness 
matrix [QJ and the damping matrix [rjJ are provided respectively by the following 

transformations, [12]: 
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[QJ= W'lQm* 


( 1 ) 


hj= 


( 2 ) 


where, [Q,] is the stiffness matrix of the composite in the material coordinate system, and 
[rj,] represents the on-axis composite damping matrix. The transformation matrix [R] and 
the matrices [Q,] and [rj,] are shown in the Appendix. The matrices have the following 

general form: 


lej - 


Q,u Qcn 0 0 0 

Q cI2 Q'22 0 0 0 ^26 

0 o Q c „ 0 0 0 | (3) 

0 0 o Q clt Q m o 

0 0 0 Q c4S Q cSS 0 

Qck Qc* 0 0 0 

and 

^cii V c j 2 0 0 0 r\ cl6 

v c2 i V c22 0 o o %26 

o o n c33 o o o (4) 

0 0 o r) c44 r) c45 0 

0 0 0 r] c54 n c55 0 

[%6i V c62 0 0 0 rj c66 \ 

Assuming viscoelastic constituents, the viscoelastic law between generalized stresses and 
strains at the ply level may be eventually expressed by the following convolute product [13], 


[V c ] = 
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a c - [QJ*e c - f‘jQ(t-T)]de(x) W 

where * indicates the convolution operation as defined above. The viscoelastic matrix 
[QJ combines both stiffness and damping matrices. The relationship is more transparent 
in the frequency domain, where the complex stiffness matrix [Q J relates the amplitudes 
of stress and strain tensors as follows: 

o c (oj) = [<?»K<«) = < w > (6) 

Hence, eqs. (5,6) represent both elastic and damping properties of the composite. The 
case of viscous damping is a special case of eqs. (5,6) 


Interlaminar Layers 

The damping polymer in each interlaminar layer is assumed isotropic, therefore each 
interlaminar layer is treated as an isotropic ply which is a trivial case of the previously 
described mechanics. Two independent elastic properties (Young’s modulus E and Shear 
Modulus G) and two independent damping properties (normal r/„ and shear r? s loss 
factors) completely characterize their viscoelastic behavior at a given frequency and 
temperature. In this manner, all independent elastic and dissipative properties of the 
composite plies and interlaminar layers are considered. 

Laminate Damping 

In order to depart from the restrictions of simplified classical laminate mechanics and 
enable the modeling of special composite laminates including the effects of thicker 
sections, interlaminar inhomogeneity, and compliant interlaminar layers, a discrete layer 
laminate damping theory (DLDT) incorporating a piecewise continuous displacement 
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field through-the-thickness is utilized [10]. Compared to other laminate theories based on 
assumed continuous displacement fields, the DLDT provided the capacity for more 
accurate strain and dissipative energy calculations in each ply, and particularly in each 
damping layer. The kinematic assumptions for the laminate theory are schematically 
shown in Fig lb. In brief, the assumed displacement field has the form suggested by 
Barbero and Reddy [14]: 

u(xyjj) = u °(xyjt) + \xy,t)F j {z) 

H 

N ( 7 > 

v(xy^) = v °(xy/)+ £ V j(xy,t)F>(z) 

H 

w(xyzj)= vv °(W) 

where, superscript o represents the uniform through-the-thickness midplane deflection, 
and u* and v* are displacements, along the x and y directions respectively, at the 
interfaces between composite plies or sublaminates (group of plies) and the interlaminar 
damping layers. F j (z) are interpolation functions. In this manner, the assumed in-plane 
displacement field is general, in that it may represent extensional, flexural, shear, and 
coupled deformations, as well as, the interlaminar shear strains through-the-thickness of 
the laminate. 

The laminate strains are directly derived 

e a (W/)= 4(w) + £ i(w) F; ( z ) 

;-i 

e ci (xy^)= e° ci (xy,t) + £ e> ci (xyf)F } /z) 
where, the midplane strains are, 


from eqs. (5) 


/-l, 2,6 
i=4,5 


( 8 ) 
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and the generalized strains are, 


e c i * 




W 



c c<5 *V 


+ V. 


( 10 ) 


The comma in the subscripts indicates differentiation. Combination of eqs. (5) or (6) and 
(8) and integration through-the-thickness of the laminate provides the following 
generalized stress-strain relations in the time domain, 


N°(f) = ( hn adz = [i(0]*e*(r)+£ [B*(f)]*«r(0 
J-w 

rf(t) = [ hl2 aF j (z)dz - [B y (01 *e°(0 " ce "(0. y= 1,..JV 
J-/V2 »-i 


( 11 ) 


where, {N°} and {A^} are the generalized stresses. The equivalent laminate relations in 
the frequency domain are: 


N 


JV*(a>) = [v4*(w)]e a (co)+£ [B**(f)]e*(o) 

m m 1 
N 

JV*( o>) = [B v (co)]e! 0 (o))+5;[I>* >, (w)]e ,,, (o)), 


j= 


( 12 ) 


where, the complex laminate stiffness matrices include both elastic and damping matrices, 


\A *] = [A] + i[A d ] 
[£*'] = [B*] * i[B’] 

[D ' im ] = [D^] + i[D?] 


In the above equation, [A d ] is the extensional laminate damping matrix, including out-of- 



plane shear terms, [B d j ] are damping coupling matrices and P d Jtn ] are flexural/shear 
damping matrices. All damping matrix expressions in eq. (13) are given in the Appendix. 
Similarly, [A] is the extensional laminate stiffness matrix, [B 5 ] are generalized coupling 
stiffness matrices, and [EH are the flexural/shear matrices given in the Appendix. 

In this context, the dissipated strain energy per unit area of the laminate Aw u which is 
defined as, 

t» L = 2ir i.TCJfeJvfe (14) 


ultimately takes the form, 


AT 


N S 


bv L • 1/2 (e e oT [A d ]€° c + 2€ c oJ Y [BfcJ + EE 

The maximum laminate strain energy per unit area w L is defined as, 


( 15 ) 


- c ^ 


(16) 


and via combination of previous equations is also expressed as a separable form of 
material properties and strains, 


N N 


w L - i/2«VK + 2efY;m< * EE <pho 

y«l ;-l 

In the above mechanics, the interlaminar compliant damping layers are considered as 


(17) 


individual plies having isotropic elastic and damping matrices [QJ and [r?J respectively. 
In this manner, the proposed laminate damping mechanics are general, as they can 
handle any composite material, number of composite plies, damping layers, and laminate 


lay-ups. 
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DYNAMIC ANALYSIS OF COMPOSITE STRUCTURES 


The developed computational mechanics for the damping analysis of composite structures 
with compliant layers, general laminations and shapes are summarized herein. The 
variational formulation of equilibrium equations for the candidate structure is as follows: 

-ffe*o c dV+f y 6u T (-pu)dV+&u l T F i + fbu?F/ls=0 ( 18 ) 

In view of the generalized laminate stress-strain relationships described in eqs. (12), the 
variational formulation in eq. (18) eventually takes the form: 

n " 

- f [&e oT N °+ 2 b^N j }dA- f { bu oT [AJu°^2bu or £ [B J m ]u J + 

JA j-l JA i-1 (19) 

If N 

£ £ bu }7 [rC\u "}<£4+8u ( r F i + fbu, T F/ls=0 

j* 1 111*1 

The first left-hand integral in eq. (19) represents both stored and dissipated strain 
energies during a virtual displacement. The second left-hand integral represents the 
kinetic energy variation in the structure, where [AJ, [B m J ] and [Dj m ] are generalized 
laminate mass matrices described in the Appendix. The third and fourth left-hand terms 
in eq. (19) represent the virtual work of discrete forces and surface tractions, respectively. 


The elastic laminate relations will result from eq. (19) if time effects are neglected in the 
first left-hand integral. The case of viscous damping is also a special case of eq. (19). In 
the frequency domain (harmonic vibration), eq. (19) results in equivalent expressions of 
complex laminate stiffness as described by eqs. (20) and (21) respectively. 

The variational relationship is discretized to an approximate equivalent system of 
dynamic equations as follows: 
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( 20 ) 


nodes 

{u w 0 } = 52 [*ji {» *» v iV*) 

i-1 
nodes 

<«*v>= E Ww 

i-1 

where, subscript i indicates the i-th node, and «, the interpolating function. The reference 
and generalized strain shape functions, [R°], and [R], respectively, interpolating the 
strains to the nodal displacements, are obtained by differentiating eq. (20) in accordance 
with eqs. (9,10). Then the damping or loss element matrix [C elj ] for nodes l and j takes 

the form: 



( 21 ) 


where A, is the area of the element. The stiffness and mass matrices of the element 
[K^j] and [M cij ] respectively, have analogous forms. 


In this context, a bilinear four-node plate element was developed. The area integration in 
eq. (21) is performed numerically. However, selective reduced integration is used for 
both shear stiffness and damping terms to avoid overstiffening of the element m the case 
of low thicknesses. Such overstiffening was observed with full integration. 


After the discretization, the system of dynamic equations has the following form m the 
frequency domain: 
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o 2 [M]{«} + i[C(<*>)]{«} + [*(w)]{w} - {*{«)} 


( 22 ) 


Viscous damping is a special case of the above equation, and results to the familiar 
dynamic system, 


[M]{u) * [£]{«} + ■ WM 


(23) 


Either one of eq. (22) or (23) may be solved directly to provide the transfer functions of 
the structure. The modal damping and natural frequencies of the structure may be 
directly obtained from the modal analysis of dynamic system (22). In a simpler form, 
however, the modal damping of the n-th mode r) n can be expressed as the ratio of the 
dissipated over the maximum modal strain energy of the structure. 


1 EAt *W. (24) 

"" ' 2« EJF*. 

where, the dissipated and maximum stored strain energies of the element for the n-th 
mode are given by: 


a*v 1/2 


I / 2 


(25) 


The frequency effects on the damping and stiffness have been presently ignored, simply 
because efficient computation may require special solution techniques which are beyond 
the scope of this paper, that is, to analyze and quantify the effect of the composite 
material parameters on the structural dynamic characteristics. However, the frequency 
effects on fibers, matrix, and interphase layers may be incorporated with minimal effort, 
and this capability of the mechanics has been already demonstrated for the case of semi 
analytical plate solutions [10]. 


12 



APPLICATIONS AND DISCUSSION 


The developed computational mechanics were applied for the analysis of damping in 
various composite plates with embedded compliant interlaminar damping layers. The 
objectives of the case studies were: (1) to evaluate the quality of the method by 
correlating predictions with both analytical and experimental results; and (2) to explore 
and quantify the potential to improve the vibroacoustic characteristics of composite 
structures by embedding compliant layers between plies of strong interlaminar 
inhomogeneity. The second group of applications demonstrates the potential to 
engineer/tailor such specialty composite laminates for superior dynamic performance. 


Materials and Assumptions 

The selected materials were 0.60 FVR T300/934 composite, and ScotchDamp ISD110 
damping film (3M Corp., St. Paul, MN) as interlaminar layers. The shear modulus and 
shear loss factor of the damping polymer were assumed 5.5 kpsi and 0.16 respectively, 
and they represent average values of the polymer at 500 Hz and 30° C as provided by 
the manufacturer’s data sheet. Due to non-existing data, the Young s modulus of the 
damping polymer was assumed 15.4 kpsi (equivalent to Poisson’s ratio of 0.3) and the 
normal loss factor was assumed to be negligible (0.016). Parametric studies have 
indicated strong insensitivity of modal damping and frequency to the dilatational 
properties of the damping polymer. Neither measured data for the damping of the 
T300/934 composite, nor for the 934 epoxy were available at the time of this report, 
therefore, the damping properties of the 934 epoxy were assumed equal to typical 
damping values of other epoxy matrices. Unless otherwise stated, the nominal thicknesses 
of the Gr/Ep plies and damping layers were 0.0058 in (0.15 mm) and 0.005 in (0.13 mm), 

respectively. 
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Verification 


[<yi/90 4 ]s Simply Supported Plate. The quality and accuracy of the developed structural 
damping mechanics were first investigated by comparisons with exact analytical solutions. 
Fig. 2 shows the modal damping and natural frequencies of a 28 cm by 28 cm (11 in by 
11 in) [0 4 /i/ 90 4 ]s simply supported plate predicted by the discrete finite element method 
and an exact plate solution previously developed based on DLDT [10], Frequency effects 
were neglected in both solutions. The finite element analysis involved a 10 by 10 uniform 
mesh subdivision. As seen in Fig. 2, there is excellent agreement between the finite 
element based and exact damping predictions. The accuracy of the predictions is slightly 
but gradually reduced at higher modes, because the mesh is too coarse for representing 
the higher mode shapes. A more refined mesh would readily improve the predictions at 
higher modes. 

To illustrate further the characteristic convergence of the element, the damping and 
natural frequencies of five mode shapes predicted with various mesh densities are plotted 
in Fig. 3. Both natural frequencies and modal damping have been normalized by the 
corresponding values of the exact solution. Clearly, both damping and frequency finite 
element predictions converge to the exact solution as the mesh is refined. Finally, as seen 
in Figs 2-3, the modal damping is typically overestimated by the discrete finite element 

model. 

Angle-ply Unsupported Plates. The predicted dynamic characteristics of two symmetric 
angle-ply "free-free-free-free" square plates with compliant damping layers are also shown 
in Figs. 4 and 5. The laminate configurations were [45 2 /-45 2 /i/45 2 /-45 2 ] s and [22.5 2 /- 
22.5 2 /i/22.5 2 /-22.5J s respectively. In-house experimentally measured damping values and 
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natural frequencies for these plates in a "unsupported" configuration [15] are also 
shown in these figures. The composite plates were fabricated in-house by hand lay-up, co- 
cured at 350° F (175° C) and 50 psi (345 KPa), and finally cut to 28 cm by 28 cm (11 in 

by 11 in). 

In both cases, the correlation between measured and predicted modal damping and 
corresponding natural frequencies of the damped plates appears very good, especially in 
the lower frequency regime. There is an overestimation of damping at higher modes 
which can be attributed to the insufficient mesh density and/or the omission of frequency 
effects. A more refined mesh and the inclusion the gradual damping reduction in ISD110 
polymer at higher frequencies will reduce the overestimation. The inclusion of frequency 
effects in the model is planned as future research. 


Predicted and measured results for composite plates of identical laminations but without 
damping layers are also shown in Figs. 4 and 5. It is pointed out, that the sole purpose of 
plotting the results for the undamped plates is for demonstrating only approximately the 
attained improvements in modal damping values when the damping layers are embedded. 
At this point the predicted and measured results for the undamped plates cannot be 
correlated, since measured damping data for the T300/934 epoxy system were not 
available at the time of this paper, and the predictions were based on assumed values. 


Damping Enhancement in Angle-Ply Laminates 


The predicted loss factors and corresponding natural frequencies of the first five modes 
of two simply-supported plates, with laminations [ 0 2 /- 0 2 /i/fV'® 2 ]s anc * 
respectively, are shown in Figs. 6-10 as functions of the ply angle 0. The modal loss 
factors and natural frequencies of plates with identical geometrical and laminate 
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configurations, but without damping layers (undamped), that is [ 0 2 /- 0 2 / 02 /- 02 ] s and [*V' 
B 2 I-Q 2 I+Q^ respectively, are also plotted in the same Figures. 

Clearly, the addition of the damping layers drastically improved the damping capacity of 
the composite plates. Most importantly, the predicted results demonstrate the potential 
for enhancing the damping of the composite structure by altering the fiber orientation 
angle. For most modes the maximum damping occurred consistently at 0=45°. Indeed in 
the considered laminate configurations, non-zero ply angles alter the anisotropy of the 
restraining sublaminates, induce additional coupling terms in their stiffness matrices and 
increase the through-the-thickness inhomogeneity of the laminate. Consequently, higher 
interlaminar shear stresses are developed in the damping layer to maintain the stress 
equilibrium, which increase the damping contributions of the damping layers. This 
damping enhancement is further signified by the observation that the predicted damping 
of the composite plates without damping layers is typically reduced for ply angles 
between 0° and 90°. Finally, the damped plates exhibit similar or higher natural 
frequencies than the undamped plates between 0 and 90 degs., hence, this enhancement 
of damping does not seem to reduce the natural frequencies and stiffness. 


As seen in Figs. 6-10, the improvements in the damping are increased at higher modes, 
for example mode 1 (1,1), mode 2 (1,2), mode 3 (1,3), and mode 5 (2,2). The numbers in 
parenthesis indicate the mode order at 0=0 degs. The uneven damping increases for 
various mode shapes indicate the interrelation between laminate configuration and mode 
shape. A significant part of the damping improvements were induced by altering the 
anisotropy in the restraining composite sublaminates. It seems though, that alternating 
fiber orientations in the plies which contain the damping layer produce some additional 
damping improvements. It is also interesting to note, that not all modes exhibited 
potential for damping improvements with off-axis plies. This is illustrated in mode 4, 
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which corresponds to mode shape (2,1) at 0°, and exhibits its higher damping at 0=0 
degs, that is, when the outer restraining sublaminates have the higher flexural stiffness in 
the direction of the flexure. 

The presented results, illustrated the accuracy and versatility of the structural damping 
mechanics in predicting the damped dynamic characteristics of specialty composite 
structures with cocured interlaminar damping layers. The applications have demonstrated 
the potential to develop "engineered" composite laminates with embedded interlaminar 
damping layers which are capable of providing higher and tailorable damping than 
traditional isotropic materials. The results finally demonstrate the existing opportunity to 
tailor such composite structures for improved dynamic characteristics in terms of 
damping and natural frequencies. 


SUMMARY 

Damping mechanics for the analysis of damping in specialty composite structures 
incorporating damping layers were developed and presented. Novel finite-element based 
mechanics incorporating a discrete layer damping theory (or layer-wise damping 
mechanics) were developed for predicting the damped characteristics of composite 
structures with embedded compliant interlaminar layers of general thickness, shape, 
lamination, and boundary conditions. 

Comparisons with exact solutions and experimental data illustrated the quality and 
accuracy of the structural damping mechanics. Case studies with interlaminar damping 
layers embedded into angle-ply composite laminates demonstrated that such angle-ply 
composite laminates and structures can combine higher and tailorable damping together 
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with improved dynamic characteristics. 
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APPENDIX 


On axis ply stiffness matrix [Q,] of a transversely isotropic composite 



Qlll Qui 

0 

0 

0 

0 ' 



Qu2 @122 

0 

0 

0 

0 



0 0 

@133 

0 

0 

0 


m= 

0 0 

0 

@144 

0 

0 



0 0 

0 

0 

@iss 

0 



0 0 

0 

0 

0 

@166 


vhere 







Q u r 



Qm - 


II 



<?»)' 

i 

f 1 

'US 

Qm‘ 

G 123 


Qm- 

G U3 

Qus “ 

°U2 


A 

= 

1-v 

U2 V t21 


Transformation matrices: 


[R] = 


m 


2 n 2 0 0 0 


n 2 m 2 0 0 0 

0 0 10 0 

0 0 0 m -n 

0 0 0 n m 


2mn 

-2mn 

0 

0 

0 


-mn mn 0 0 0 m 2 -n 


m = cos0, n = sin0 


[R]-‘ = [R(-e)] 

On-axis damping matrix: 


( 26 ) 


(27) 


(28) 


(29) 
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[*7/] - diag { rjjj, i) l2 , rj l4 , t) 15 , r\ l6 } 
Generalized damping matrices 


( 30 ) 


KJ= 2^ E [*“ 

t-i J*»-i 


27T £ f Ak 
*-l J *k-1 

*1 . 

z> =1,2,6 

2t e p mMnxwkyk 

k - 1 J "k-1 
*1 . 

i/i =4,5 
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Generalized stiffness matrices 
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Generalized mass matrices 
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*-i 

"i . 

pji- E [ [*» wy^-CO* 

Jk-l JA k-’ 
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(b) Kinematic assumptions of discrete layer laminate theory. 


Figure 1 Specialty laminates with interlaminar damping layers, (a) Typical laminate 
configuration; (b) Kinematic assumptions of discrete layer laminate theory. 
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Figure 2 Modal characteristics of a (0 4 /i/90 4 )s T300/934 simply supported plate with 
ISD110 damping layers. Numbers in parenthesis indicate mode order of 
exact solution. 
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Predicted and measured modal characteristics of a [45 2 /-45 2 >^/45 2 /-45 2 ] s 
T300/934 simply supported plate with ISD110 damping layers. 
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Figure 5 Predicted and measured modal characteristics of a [22.5 2 /-22.5 2 /i/22.5 2 /- 
22.5 2 ] s T300/934 simply supported plate with ISD110 damping layers. 
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THETA, DEGS 


Modal loss factor and natural frequency of angle-ply T300/934 simply 
supported plates; mode 1. 
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Modal loss factor and natural frequency of angle-ply T300/934 simply 
supported plates; mode 4. 
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THETA, DEGS 



Figure 10 Modal loss factor and natural frequency of angle-ply T300/934 simply 
supported plates; mode 5. 
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